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Abstract. The cluster Abell 2104 is one of the lowest red- 
shift clusters [z — 0.153) known to have a gravitational 
lensing arc. We present detailed analysis of the cluster 
properties such as the gravitational potential using the 
X-ray data from ROSAT (HRI) and ASCA, as well as 
optical imaging and spectroscopic data from the CFHT. 
The cluster is highly luminous in the X-ray with a bolo- 
metric luminosity of L x ~ 3 x 10 45 ergs s _1 and a high gas 
temperature of ~ 10.4 keV. The X-ray emission extend- 
ing out to at least a radius of 1.46 Mpc, displays signif- 
icant substructure. The total mass deduced from the X- 
ray data under the assumption of hydrostatic equilibrium 
and isothermal gas, is found to be M tot (r < 1.46Mpc) ~ 
(8.0 ± 0.8) x lO 14 Af . The gas fraction within a radius of 
1.46 Mpc is ~ 5 — 10%. The cluster galaxy velocity distri- 
bution has a dispersion of 1200±200 kms -1 with no obvi- 
ous evidence for substructure. The total mass within 1.46 
Mpc, deduced from Jean's equation using the observed 
galaxy number density distribution and velocity disper- 
sion, is found to be ~ 6.8 x 10 14 Af Q to - 2.6 x 1O 15 M 
marginally consistent with the X-ray deduced total mass. 
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1. Introduction 

Clusters of galaxies are the largest bound systems in the 
Universe, and as such they are the largest objects where 
detailed studies of their gravitational potential are pos- 
sible. Given their large sizes, 3 to 6 Mpc in extent, they 
are also thought to be representative of the Universe in 
terms of the baryonic fraction which is directly related to 
the density of the universe and the predictions of the Big 
Bang nucleosynthesis theory. Studies so far have found 
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that the baryonic fractions in clusters favour a low mat- 
ter density universe given the predictions of baryon den- 
sities given by the nucleosynthesis theory (e.g. White et 
al. 1993). Recently, detailed and independent estimates of 
cluster total mass distributions have become available; the 
mass-tracers used and the observational techniques em- 
ployed can be summarised as follows: 

— Cluster Galaxies : these have a long tradition of provid- 
ing mass estimates via application of the Virial The- 
orem to the observed dispersion in their radial veloci- 
ties. The method rests upon the assumption that the 
galaxies are in dynamical equilibrium. 

— Hot Intracluster Gas : as well as being an important 
mass component of clusters, its X-ray emission pro- 
vides an ideal tracer - through the hydrostatic equa- 
tion - of the total underlying mass. The assumption 
that the gas is in hydrostatic equilibrium with the clus- 
ter's gravitational potential is thought to be reason- 
ably secure for the central few Mpc (Evrard et al. 1996 
and Schindler 1996) and the gas density and tempera- 
ture profiles required to solve the hydrostatic equation 
are readily available from the X-ray data. 

— Gravitational Lensing : here the lensing action of the 
cluster on background sources, as revealed in deep high 
resolution imagery (Tyson et al. 1990; Fort & Mellier 
1994 and references there in), is used to provide a di- 
rect measure of the shape and depth of the cluster 
potential and hence the projected mass distribution 
(Kaiser & Squires 1993, Broadhurst et al. 1995 etc.). 
Unlike the first 2 methods, this approach is not reliant 
upon assumptions of hydrostatic or dynamical equilib- 



For detailed studies in the X-ray and optical, we need 
a nearby cluster, though gravitational lensing effects are 
diminished for low redshift clusters. An ideal cluster for 
this kind of detailed and independent estimates of mass 
distributions, would be one of the lowest redshift clusters 
with obvious lensing effects such as a giant arc. In this 
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paper, we will analyse the X-ray and optical data for one 
of the nearby lensing clusters. 

Abell 2104 is a rich cluster (richness class 2) at a red- 
shift of 0.155 (Allen et al. 1992). It was found to have 
a high X-ray luminosity from the ROSAT all-sky survey 
data (Pierre et al. 1994). Subsequent optical followup ob- 
servations with the CFHT revealed an arc embedded in 
the halo of the central cD galaxy 7'.'2 away from the cen- 
tre (Pierre et al. 1994). The arc spans 10" in length and it 
is amongst the reddest known arcs. Fig. || shows a close up 
picture of the arc. Given the small arc radius, it is impor- 
tant to have a high resolution X-ray observation with an 
instrument such as the ROSAT/HRl to probe the gravi- 
tational potential within the arc radius. 

The optical data including photometry and spec- 
troscopy will be analysed in Sec. 2. The spatial and spec- 
troscopic analysis of the X-ray data from ROSAT and 
ASCA will be given in Sec. 3. The independent mass es- 
timates using different methods as well as a comparisons 
will be given Sec. 4. 

Throughout the paper we adopt a cosmological model 
with H a = 50 kms _1 Mpc~ 1 , Sl = 1 and A = 0. Celestial 
coordinates are in J2000. 
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Fig. 2. A close up image of the central regions of 
Abell 2104 showing the giant arc. 



2. Optical data 

2.1. Observations 

The data were collected in 4 nights at the 3.6 m CFHT 
Telescope in May 1993. Two 10 minutes exposures in B 
band and two 15 minutes exposures in R band were ob- 
tained. Exposures of 30 to 55 minutes per spectroscopic 
mask was obtained for 3 separate masks, each containing 
about 30 slits (Fig. |). The focal reducer MOS/SIS to- 
gether with CCD Lick2 (2048 x 2048 pixels of 15 fj,m) were 
used during the run. This CCD is a thick device having a 
quantum efficiency of ~ 10% in the blue. The observing 
configuration provides a pixel size of 0.314" over a field of 
view of about 10' x 10'. The overall image quality was good 
(stellar FWHM ~ 0.9") although some optical distortions 
were conspicuous near the edges of the images due to the 
optics of the focal reducer. 

2.2. Photometric analysis 

The B and R frames were prepared using standard pre- 
reduction techniques. Since there were only 2 frames per 
filter, cosmic rays were removed by taking the lower pixel 
value in cases where a pixel in one frame is significantly 
higher than the corresponding pixel in the other frame. 
The photometric analysis was performed by means of the 
SExtractor package (Bertin & Arnouts, 1995) in the same 
way as Pierre et al. (1997), but adapted to our data. The 
images were first slightly smoothed to give the same PSF 
in B and R frames, then the background was estimated us- 
ing a 64 x 64 pixel mesh. Source detections were claimed 



if at least 9 adjacent pixels were above a threshold cor- 
responding to 1.5 times the local noise level. The CCD 
Sequence in M 92 (Christian et al., 1985) observed dur- 
ing the same run was used for photometric calibration. 
Stars VCS1, A, B (probably variable) had to be removed 
because of obvious inconsistencies. Estimates of the pho- 
tometric errors were taken directly from the SExtractor 
analysis, and are less than ~ 0.1 for R< 22.5 and less 
than - 0.2 for B< 23.5. 

The catalogue is estimated to be complete to R = 22.5 
and B = 23.5. On inspection of the detected objects above 
the completeness limit, we found those objects with a SEx- 
tractor classification < 0.15 may be assumed to be galax- 
ies, i.e., 275 objects. Changing the threshold does not af- 
fect the outcome significantly because most of the galaxies 
are well separated from stars (3/4 of the objects fall below 
0.05 or above 0.95). 

Fig. H shows the colour magnitude diagram for all the 
galaxies detected in the R-frame, and the corresponding 
magnitudes in the R and B bands were measured within 
the same apertures. The band of E/S0 sequence galaxies 
is discernable in Fig. ||; the spectroscopically confirmed 
cluster members are shown to fall mostly on the E/S0 
sequence confirming that a large fraction of the galaxies 
on the E/S0 sequence belongs to the cluster. The mean 
error in B-R colour is < 0.08. 

2.3. Spectroscopy 

Grism O300 was used for the spectroscopy. It has a zero 
deviation at 5900 A, covers approximately 4700-7900 A, 
and gives a dispersion of 3.59 A/pixel (0.314"). The slit 
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Right Ascension (JS000) 

Fig. 1. Optical field of Abell 2104, observed at the CFHT in R band. Overlaid are the ROSATHRJ contours with levels 
(1.7, 2.0, 2.6, 3.2, 3.8, 4.3, 4.6) x 10~ 6 counts s _1 arcsec~ 2 . The X-ray image was rebinned into 2 pixels and smoothed 
with a 10 Gaussian. 



has a width of 2", i.e. 6.4 pixels, yielding a resolution of 
~23 A FWHM. Since there was only 1 frame per mask, 
cosmic rays were picked out individually by eye and re- 
placed by the median of the surrounding pixels. The in- 
ternal Helium and Argon lamps was used for wavelength 
calibration. The subsequent reduction was performed as 
described in Pierre et al. (1997). Redshifts were measured 
by a cross-correlation method implemented in the MIDAS 
environment following Tonry and Davis (1979). The cross- 
correlation results for each spectrum were checked inde- 
pendently by eye. 

The results from the cross-correlation analysis for all 
spectra are presented in Table [j[ Heliocentric correction 



has not been applied, but is negligible at this resolu- 
tion. The absolute error in the velocity calibration is 
- 200 kms" 1 . 

As a first guess, galaxies are considered to be cluster 
members if they lie within 3000 kms -1 of the central cD 
galaxy which selects 47 (the main sample) out of the 60 
galaxies. This procedure eliminates most of the foreground 
and background galaxies without affecting the dispersion 
measurements significantly. If we relax the velocity con- 
straint and apply the usual 3a— clipping technique then we 
have 51 cluster members (the extended sample). The clus- 
ter redshift distribution for both samples is displayed in 
Fig. ||. The histogram includes all galaxies in the redshift 
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Fig. 3. Finding chart for galaxies with measured redshifts. The reference numbers are the same as in Table |[ The 
non-member galaxies are marked with a bracket around their reference number. 
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Fig. 4. A colour magnitude diagram for all galaxies de- 
tected in the R frame. The filled circles are for spectro- 
scopically confirmed cluster members, and the crosses are 
for the non-members. The solid line gives the completeness 
limit. 



range z ~ 0.135 - 0.175 in Table g and a Gaussian corre- 
sponding to the velocity distribution of the main sample. 
The bi-weighted mean and scale for the main sample are 
z = 0.1532t8;gg{$ and a = 1148±^§° kms^ 1 correspond- 
ingly; and z = 0.1538±g;ggg| and a = UOll^kms- 1 for 
the extended sample. It is difficult to find an objective 
criterion for deciding which galaxies are cluster members. 
Even with the sophisticated weighting scheme employed 
by Carlberg et al. (1997), the determination of the weight 
for each galaxy is still subjective. In Table g, we have 
marked only the galaxies from the main sample as cluster 
members. 

For the main sample we have enough redshifts to test 
whether or not the galaxy velocities are drawn from a 
Gaussian distribution applying various statistical tests for 
normality (e.g. D'Agostino & Stephens 1986; ROSTAT 
- Beers et al, 1990; Bird & Beers, 1993). As a result 
Anderson-Darling test (A 2 ) accepts the hypothesis for nor- 
mal distribution at 90% significance, the combined skew- 
ness and kurtosis test (Bl & B2 omnibus test) at 97% 
level and the alternative shape estimators, asymmetry in- 
dex and tail index based on order statistics, were found to 
be —0.21 and 0.98 respectively which also show that the 
velocity distribution is drawn from a Gaussian. 

We can obtain a conservative estimate of the errors 
on the velocity dispersion by comparing the dispersion 
from the extended and main samples. When we take into 
account of the uncertainties in cluster membership, a more 
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Fig. 5. Cluster galaxy redshift histogram (bin size Az = 
0.0025). The galaxies considered to be cluster members are 
marked as solid histogram while the dotted histogram are 
for the extended sample. The dashed curve is a Gaussian 
with parameters corresponding to the velocity distribution 
of the galaxies in the main sample. 



conservative estimate of the errors should give the velocity \ ^ | 

dispersion as 1200 ± 200 kms -1 . 

We also investigated the presence of substructures in 
(a, <^ z) space but no obvious signal was detected (see Fig. 6. Wedge diagrams showing the distribution of galax- 
Fig. g). More redshifts are required for a proper statistical j es w ith measured redshifts 
analysis. 



3. X-ray Data 

We have observed the cluster with the ROSAT HRI and 
the ASCA GIS and SIS detectors. The HRI has a high 
spatial resolution of ~ 5", which provides a high resolu- 
tion X-ray surface brightness profile, but it has no energy 
resolution. ASCA on the other hand has a low spatial res- 
olution (~ 3') but relatively high energy resolution and 
high sensitivity in the energy range 1-10 keV, which pro- 
vides a reliable gas temperature measurement for clusters 
of galaxies. 

3.1. Spectral analysis 

The cluster was observed with ASCA using both detec- 
tors of the Gas Scintillation Imaging Spectrometers (GIS) 
and Solid-state Imaging Spectrometers (SIS) in February 
1996. The SIS detectors were operated in 1-CCD mode. 
The data was screened and cleaned according to the stan- 
dard procedures recommended (The ABC guide to ASCA 



data reduction) . The spectra were extracted from the cen- 
tral ~ 6.5' radius from the GIS2 and GIS3 detectors, ex- 
cluding one discrete source. Similarly, spectra were ex- 
tracted from the central ~ 3' radius from the SIS0 and 
SIS1 detectors. A standard blank-sky exposure screened 
and cleaned in the same way as the cluster field was used 
for background subtraction by extracting a background 
spectra from the same region on the detector as the clus- 
ter spectra. The spectra were grouped into energy bins 
such that the minimum number of counts before back- 
ground subtraction was above 40, which ensures that x 2 
statistics would still be valid. The 4 spectra from each de- 
tector were simultaneously fitted with a Raymond-Smith 
thermal spectra (Raymond & Smith 1977) with photoelec- 
tric absorption (Morrison & McCammon 1983) from the 
XSPEC package (Fig. ^) . We adopted the abundance table 
with the relative abundance of the various elements from 
Feldman (1992). The free parameters were the gas temper- 
ature (Tg), Galactic neutral hydrogen absorption column 
density (N(H)), metal abundance (abund) and the emis- 
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Fig. 7. ASCA spectra from the 4 detectors GIS2, GIS3, 
SISO, SIS1. The solid curves show the simultaneous fit to 
all 4 spectra using a Raymond-Smith model with photo- 
electric absorption. The model fits also take into account 
of the instrumental responses of the individual detectors. 

sion integral. All 4 spectra were to have the same value 
for the free parameters except for the emission integral, 
since the GIS and SIS PSF were different and the extrac- 
tion regions were smaller for the SIS spectra compared to 
that of the GIS. The two GIS spectra were assumed to 
have the same emission integral but different from the SIS 
emission integrals. Results of the best simultaneous fit to 
the 4 spectra along with fits to the individual spectra are 
tabulated in Table [l]. Only data in the energy range where 
the effective area of the detectors arc > 10 cm 2 were used 
for the spectral fitting, i.e. 0.6-7.5 keV for SIS data and 
0.85-10.0 keV for GIS data. 

The neutral hydrogen column density derived from 
the ASCA data were 2 times larger than the N(H) (= 
9.25 x 10 20 cm 2 ) measured from radio data by Starck 
(1992). If we try to fix N(H) to the value determined by 
Starck (1992), then there is obvious discrepancy between 
the model spectrum and the SIS data below 1 keV. Un- 
fortunately, there is no PSPC data available for this clus- 
ter to place definitive constraints on the N(H) value. It 
is possible that there is a local over-density of absorbing- 
neutral gas along the line-of-sight to the cluster, though it 
is more likely to be a calibration error for the SIS detector. 
Calibration of the low-energy part of the SIS detector is 
known to produce erroneous results such that it favours a 
high N(H) inconsistent with PSPC results (Schindlcr ct al. 
1998 & Liang et al. 2000). In view of the possible calibra- 
tion error for the SIS, the data were also fitted with the 
above models with a fixed N(H) given by Starck (1992) 
by excluding the SIS data below 1 keV. The temperature 
thus deduced was significantly higher than before. In the 
following studies, we will adopt these parameters deduced 
from a simultaneous fit of data from the GIS detectors 
in the energy range 0.85 to 10 keV and the SIS detectors 
between 1 and 7.5 keV. 



3.2. ROSAT HRI data 

The cluster was observed by the ROSAT HRI in February 
(7.6ksec) and August (36ksec) 1996. The X-ray centroid 
was found to be 15:40:08.1 -03:18:17, which is - 1" from 
the position of the cD galaxy 15:40:07.96 -03:18:16.7. 
The positional error for the X-ray centroid is ~ 5", hence 
the small apparent displacement between the cD position 
and the X-ray centroid is insignificant. The X-ray surface 
brightness was obtained by extracting the photons in a ra- 
dius of 7' and the background was extracted from an annu- 
lus of 8' — 10' radius from the X-ray peak. Discrete X-ray 
sources were excluded from the extraction. There were 8 
discrete X-ray sources in the HRI image. Fig. [l] shows the 
X-ray contours overlaid on the optical image of the clus- 
ter field. The X-ray image show significant substructure in 
the centre with an overall elliptical appearance. The dis- 
crete X-ray sources at 15:40:07.2 —03:19:53 is embedded 
in the cluster emission. The relative astrometry between 
X-ray and optical was checked using 4 of the discrete X- 
ray sources that had clear optical identification. The X-ray 
positions had a maximum displacement of ~ 2" relative to 
the optical coordinates. The X-ray contours were adjusted 
to the optical coordinate system using the 4 discrete X- 
ray sources, which gave a relative astrometric accuracy of 
~ 0.5" between optical and X-ray coordinates. 

A radial average of the X-ray surface brightness for the 
cluster is shown in Fig. ||. A best fit f3 profile (Cavaliere 
and Fusco-Femiano 1976) 

S x (r)^S [l + (-f]- 3P+1/2 (1) 

convolved with the instrument PSF is shown as a solid 
curve superimposed on the data. The best fit gave /? = 
O.SOIq'q 2 and 9 — 51^|". The uncertainties quoted are la. 
The total X-ray luminosity within the central 7' radius is 
L x ~ 8.95 x 10 44 ergs s" 1 in the ROSAT band of 0.1- 
2.4 keV, assuming N(H) = 9.25 x 10 20 cm 2 , kT g = 10.4 
keV (or 1.2 x 10 8 K), and an abundance of 0.22. The X- 
ray luminosity thus deduced is consistent with that es- 
timated from the ROSAT all-sky survey (Pierre et al. 
1997). The corresponding bolometric X-ray luminosity is 
L x ~ 3 x 10 45 ergs s _1 . The central electron density was 
thus derived to be n e fl ~ 5.92 x 10~ 3 cm -3 . The central 
cooling time for this cluster is t coo i ~ 10 10 yr, greater than 
a Hubble time. 



4. Analysis 

While the X-ray image show significant substructure in 
the cluster indicating deviations from hydrostatic equilib- 
rium, the cluster total mass deduced from assumptions 
of dynamical equilibrium are still reliable, as is shown 
by numerical simulations (Evrard et al 1996 & Schindler 
1996). Under the assumption of hydrostatic equilibrium 
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Table 1. Results on spectral fit to ASCA data 
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<-°°-0.52 




11.04+^ 


10 S6 +U b4 


abund 


21+° 11 


22+° 11 


0.32 ±0.08 


0.20 ±0.13 


0.23 ± 16 


0.22 ± 0.07 


N(H) 

x 2 


15.3 ±4.7 


22.5 ± 2.2 


20.8 ± 1.6 


9.25 


9.25 


9.25 


0.55 


0.75 


0.91 


0.56 


0.90 


0.95 



Notes: 

kT 9 - the gas temperature in keV; 

abund - the fractional solar metal abundance; 

N(H) - the neutral hydrogen column density in units of 10 20 cm 2 ; 

X 2 - reduced \ 2 ■ 

col. 2 - fit to the combined GIS data; 
col. 3 - fit to the combined SIS data; 

col. 4 - simultaneous fit to GIS2, GIS3, SIS0 & SIS1 spectra; 

col. 5,6,7 - same as col. 2,3,4 respectively, but N(H) was fixed to the radio value and SIS data below 1 keV were not used. 
The quoted errors for each parameter correspond to the 90% confidence range. 
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Fig. 8. X-ray surface brightness from HRI data. The solid 
curve gives the best (3 profile after convolution with the 
HRI PSF. The dashed horizontal line segment indicates 
the background level. 

and spherical symmetry, the cluster total mass is directly 
related to the intracluster gas properties as: 

rkT g (r) d In n e (r) d\uT g (r) 



M tot {r) 



-(- 



dim 



din? 



) 



(2) 



[im v G 

In general, a good fit can be found for the X-ray surface 
brightness distribution using the parametrisation given in 
Eqn. |], which in turn gives the gas density as follows if 
the gas is isothermal: 



^ = [l±(r/r ) 2 r 3/3/2 



n e ,o 



(•3) 



Hence, the gravitational potential is given by 



(r) - 0o 



3a 2 



2 



(4) 



where cr 2 = (3kT g / /im p and the total mass is given by 
M tot {r) = (^)- (r/r ° )3 



G 'l + (r/r Q y 



(5) 



The lensing effects of the background galaxies by the 
cluster gravitational field is directly related to the 2-D pro- 
jection of the total mass density. In this case, the projected 
total mass density is given by 



V 2D M _ 3^ 2±(r/r )^ 



4Gr [l±(r/r ) 2 ] 3 / 2 ' 



(6) 



If we consider galaxies as test particles in the clus- 
ter potential well, then Jean's equation for a collisionless, 
steady state, non-rotating spherically symmetric system 
gives 



M tot {r) 



ra 2 r (r) d\nn gal (r) dinger) 



G 



d\n 1 



where n ga i is the spatial galaxy number density, (3t is 
the anisotropy index and oy is the radial velocity dis- 
persion. The spatial galaxy number density is related to 
the observed 2-D projection of the galaxy number density 
through the Abel inversion given by 



n gai(r) 



1 

27T 



dtf°(R) d R 



dR 



(8) 



cr,. and Pt are related to the observed line-of-sight velocity 
dispersion cr; through 



f°° f? 2 
S^(i?)cr 2 ( J R) = 2 / n gal (r)a 2 r (r)[l-\(3t} 



rdr 



Vr 2 - R 2 



(9) 
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In the simple case, where the galaxy orbits are 
isotropic, Eqn. is equivalent to Eqn. || with of replaced 
by kTg/nm p . 

If we make a further simplification by assuming that 
not only the gas but also the galaxies are isothermal, i.e. 
oy(r) is a constant, then we have 



n e {r) = , n ga i(r) 0a 



n e ,o 



l-gal.O 



(10) 



where (3 8 = fj,m p a 2 /kT g . Given the above parametrisation 
for the X-ray surface brightness and the resultant expres- 
sion for n e given by Eqn. |3|, we deduce the spatial galaxy 
density distribution as 



rigalfi 



= [l + (r/r ) 2 



(11) 



where a = 3/3/ 2(3 s . The observed line-of-sight velocity dis- 
persion is trivially given by a bs = o~ r and a = Z0 2 . /2a 2 bs . 

Alternatively, if we simplify the case by assuming that 
the galaxy density distribution follows that of the total 
mass, i.e. mass-follows-light, then from Jean's equation 
(Eqn. |7|) we see that the galaxies can not be isothermal if 
the gas is isothermal and the X-ray surface brightness is 
parametrised as in Eqn [j]. The radial velocity dispersion 
is given by 



a 2 (r) 



1 + (r/r ) 2 /2 
l + (r/r ) 2 /3 ( 



(12) 



where again a 2 = /3kT g / /j,m p . The line-of-sight velocity 
dispersion 07 can be deduced from Eqn. [5| However, the 
measured velocity dispersion is an average of c; within a 
certain radius: 



(< R) = i [ cw 



(13) 



In the case of Abell 2104, we have the observables 
S X (R), T g , ?? g ai{R), (Jobs- Since the ASCA PSF was too 
poor to deduce a meaningful temperature profile, we will 
assume that the gas is isothermal for the time being. In the 
following section we will study the cluster total mass de- 
duced from the various methods and examine their consis- 
tency using the simple parametrised /3-model given above. 

Mass estimate from optical data 

The projected galaxy density distribution is consistent 
with a wide range of models. The following family of 
parametrised functions 



^gal(r) = £ S ai,o[l + (r / V ) 2 



(14) 



were fitted to the projected galaxy density distribution af- 
ter background subtraction using the density of galaxies in 
the annulus 220" to 240" as background. If we fix the core 
radius to the X-ray determined value of 51", then we found 



; - v 



■ 20 - \ \ 




50 100 150 200 250 

radius (arcsec) 

Fig. 9. The radially averaged galaxy number density dis- 
tribution (with background) . The curves show a number of 
statistically consistent model fits to the data points. The 
short-dashed, solid and long-dashed curves corresponds to 
the k = 1/2, 1,3/2 cases of the family of curves given by 
Eqn. The dotted curve shows the 2D projection of the 
Navarro model (Navarro et al. 1996). 



the best fit to be k — 1 (x 2 — 9.3 with 10 degrees of free- 
dom), though k = 1/2 to 3/2 were also statistically con- 
sistent with the observed data. Note that k = 1/2,1,3/2 
corresponds to spatial galaxy distributions of the form 

1,3/2,2 respectively. The pro- 



given in Eqn. 11 with a 

jected total mass density distribution given by Eqn. ^|was 
also statistically consistent with the projected galaxy den- 
sity distribution (x 2 = 9.66 with 10 degrees of freedom), 
which means mass-follows-light is not excluded. The 2D 
projection of the functional form (r/r0) _1 [l + (r/r0)]~ 2 
(Navarro et al. 1996) was also found to be statistically 
consistent with the observed galaxy distribution. The pro- 
jected galaxy distribution is shown in Fig. ^| along with the 
various model fits. The observed galaxy density distribu- 
tion is still declining towards the edge of the image indi- 
cating a wider field is needed to reach the true "edge" of 
the cluster. The X-ray data show that the cluster extends 
at least out to a radius of 7' which is beyond the optical 
field of view for the current observation. A wider field of 
view would help to reject some of the above models. 

If we estimate the total mass distribution from the 
galaxy density distribution and velocity dispersion assum- 
ing that the galaxies are isothermal, then a 2 = {2a/2>)a 2 bs 
where the observed data give a ~ 1 to ~ 2 and <r b s = 
1200 ± 200 kms" 1 , implying that er ~ 823 to - 1625 
kms" 1 . Thus from Eqn. ^ the total mass is between 3.5 x 
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1O 14 M and 13.4 x 1O 14 M within a radius of 220" (or 
0.76 Mpc), and between 6.8 x lO 14 Af and 2.6 x 1O 15 M 
extrapolating to 7' (or 1.46 Mpc). Note that optical data 
alone does not constrain the mass very well, even under 
assumptions such as isothermality of the galaxy distribu- 
tion and isotropy of the orbits. 

On the other hand, if the galaxy distribution is not 
isothermal but follows that of the mass then the measured 
velocity dispersion implies that gq ~ 1010 ± 165 kms -1 
from Eqn. |l| and the total mass is - (10.2±3.7) x 1O 14 M 
within a radius of 7' (or 1.46 Mpc). 

4-2. Mass estimate from X-ray data 

The values of (3 ~ 0.5±g;g|, r ~ §l±l" and T g ~ 10.4±0.6 
keV have been determined from spatial analysis of the HRI 
data and the spectro- analysis of the ASCA data respec- 
tively. Thus from the X-ray data, Cq = 0kT g / fim p implies 
a (To ~ 895 ± 45 kms -1 and a X-ray deduced total mass 
of - (8.0 ± 0.8) x 1O 14 M out to a radius of 7' (or 1.46 
Mpc). 

Note that if the galaxies are isothermal, then the X- 
ray deduced mass is consistent with the optically deduced 
mass (or generalised "Virial" mass) if a ~ 1. The X-ray 
and optical data are also marginally consistent if mass- 
follows-light. 

The total gas mass within 7' was found to be ~ 
7.8 x 1O 13 M which gives a gas fraction of ~ 10% com- 
pared to the X-ray deduced mass, but 5-10% compared 
to the dynamically deduced mass. The gas fraction within 
a radius of rsoo = 1.14 Mpc (where the over-density is 500 
times the critical density of the Universe) is ~ 8%, which 
is lower than the average gas fraction of (20 ± 1.9)% for 
nearby hot (kT g > 4keV) non-cooling flow clusters (Ar- 
naud & Evrard 1999). The gas fraction within a radius 
of 1.46 Mpc gives a lower limit to the baryonic fraction. 
Since the baryonic matter density predicted from the Big 
Bang nucleosynthesis gives fib ~ 0.04 — 0.06 (Walker et 
al. 1991) from the measured light element abundance, the 
lower limit of the baryonic fraction of this cluster is thus 
consistent with Qm < 1. 

5. Discussions 

For the above simple models, we have shown that the X- 
ray deduced mass is consistent with that from the optical 
data over the scale of 1-3 Mpc under the assumptions 
of dynamic equilibrium. In a recent paper by Lewis et 
al. (1999), they also found the X-ray and dynamically de- 
duced mass were consistent for a sample of CNOC clusters 
at z ~ 0.3. 

On the other hand, in a study of a sample of clus- 
ters with giant arcs, Allen (1998) found that the X-ray 
deduced mass was consistent with the position of the gi- 
ant arcs for cooling flow clusters but 2 — 3 times smaller 
than the lensing mass for non-cooling flow clusters. This 



was then explained as a direct consequence of the the- 
ory that cooling flow clusters were dynamically more re- 
laxed than non-cooling flow clusters since cluster mergers 
would certainly disrupt a cooling flow. The cooling time 
for Abell 2104 is t coo u ng ~ 10 10 yr at the centre, thus there 
is no evidence for a cooling flow in this cluster. Pierre et 
al. (1994) found a red tangential arc 7'.'2 from the centre 
of the cD galaxy (see Fig. |^). They found that the pro- 
jected mass within the arc to be 6 x lO 12 Af . Here we 
examine if the arc feature is consistent with the simple 
cluster potential deduced from the X-ray data. Since the 
projected density must reach the critical value at 7'.'2, it 
requires oo ~ 1380— 1175 kms -1 for an arc redshift in the 
range z arc ~ 0.5 — 3 assuming the potential is spherically 
symmetric. However, the X-ray data gave oo ~ 895 ± 45 
km s" 1 apparently inconsistent with the lensing deduced 
value, indicating that in this very simplistic model the X- 
ray mass within the arc radius appears to be ~ 1.5 — 2 
times smaller than needed to produce the giant arc. Our 
result appears to be consistent with the results of Allen 
(1998). However, since the model we have adopted so far 
is very simple and the arc radius is relatively small (7'.'2), 
it is premature at this stage to suggest that the lensing 
results are inconsistent with the X-ray data under the as- 
sumptions of hydrostatic equilibrium and isothermal gas. 
As it was pointed out in Pierre et al. (1994), the small 
arc radius is an indication that the local cD potential is 
probably as important as the global cluster potential in 
forming the arc feature. Indeed for most clusters with gi- 
ant arcs, the arc radii are barely larger than the PSPC 
resolution and probably a few times larger than the HRI 
resolution, hence an inconsistency between X-ray deduced 
mass from simple models and that of the strong lensing 
deduced mass are not sufficient to prove that the cluster 
is not in dynamic equilibrium. An alternative explanation 
for the results of Allen (1998) could be that the cooling 
flow clusters are well modelled by a cluster potential simi- 
lar to the type given by Eqn. |J, but non-cooling flow clus- 
ters have a different shape of gravitational potential, e.g. 
a mass profile that has a broad component in the outer 
parts of the cluster (e.g. Gioia et al. 1998). It would be 
difficult for the HRI to reject a model of this kind since 
it has a high background level and it would be easy to 
"hide" faint diffuse emission at large radii. In our study of 
Abell 2104, the current optical image does not extend to 
the extent of the X-ray emission, thus we need wide-field 
imaging to find out the true extent of the cluster. 

The X-ray emission in the centre of the cluster shows 
strong ellipticity, the effect such asphericity has on the 
mass estimates needs to be addressed since the mass es- 
timates given above were calculated under the assump- 
tion of spherical symmetry. Neumann & Bohringer (1997), 
estimated the effects of asphericity on mass estimates 
of CL0016+16, and found that the total mass was only 
changed by ~ 2% when the ellipticity was taken into ac- 
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count. The ellipticity demonstrated in the X-ray image of 
Abell 2104 is no stronger than that of CL0016+16. 

So far we have only considered the isothermal gas mod- 
els, but the total mass given by Eqn. |^ is more sensitive 
to T g than n e . It is necessary to explore models with a 
temperature gradient. Markevitch et al. (1998) found an 
almost universal decrease in temperature in the outer re- 
gions over a radius of 0.3 to 1.8 Mpc in a sample of 30 
nearby clusters (0.04 < z < 0.09). They found that for 
a typical 7 keV cluster, the observed temperature profile 
can be approximated by a polytropic equation of state 
with 7 ~ 1.2 — 1.3. If we assume that Abell 2104 has a 
similar large scale temperature profile, then we can quan- 
tify the mass ratio between the polytropic and isothermal 
models as 



Ml^jr) n e (r) 



7-1 



(15) 



Since the X-ray emissivity has only a weak dependence on 
T g over the 1-10 keV range (only a 10% change), the X-ray 
surface brightness varies insignificantly with T g . We can 
then safely take the gas distribution as determined from 
the isothermal case (i.e. Eqn. |). Thus at 7' radius (1.46 
Mpc), a model with such a temperature gradient would 
give a mass that is ~ 0.6 times smaller than the isother- 
mal case. This would cause the X-ray deduced mass to be 
strongly inconsistent with the dynamically deduced mass 
unless of increases with radius in a similar manner as T g . 
Note that a temperature profile that decreases with the 
radius would also increase the total mass in the inner clus- 
ter regions compared to the isothermal model, and thus 
alleviate the discrepancy between the X-ray mass within 
the arc radius and the position of the giant arc. Fig. [To] 
shows the range of mass profiles deduced from the various 
methods and models discussed in the paper. 

6. Conclusions and Future prospects 

The rich cluster Abell 2104 at a redshift of z = 0.1533 
was found to have a high X-ray luminosity (~ 9.0 x 10 44 
ergs s -1 in [0.1-2.4] keV) and temperature (10.4±0.6 keV) 
from ROSAT HRI and ASCA data. The central cooling 
time, t coo i ~ 10 10 yr for this cluster indicates the absence 
of a cooling flow. The galaxy velocity distribution showed 
that the cD galaxy was at rest at the bottom of the cluster 
potential. The X-ray image shows significant substructure 
in the centre of the cluster and an overall elliptical ap- 
pearance. It appears that the cluster has not yet reached 
dynamical equilibrium. 

As shown in Evrard et al (1996) and Schindler et al. 
(1996), the total mass deduced from assumptions of dy- 
namical equilibrium are not significantly different from 
the true values. The total mass deduced from the X-ray 
data assuming hydrostatic equilibrium is consistent with 
the dynamic mass deduced from Jean's equation. How- 
ever, the current data on the projected galaxy density 
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0.01 



L. 

0.1 
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Fig. 10. A comparison of 3D total mass distribution de- 
rived by the various methods. The solid curves give the 
range of X-ray deduced mass for isothermal gas; the dot- 
ted curves give the range of mass for the polytropic model. 
The dashed curves give the range of dynamic mass derived 
from the galaxy density distribution and velocity disper- 
sion. The curves are plotted only for regions were data is 
available. The two stars show the range of mass estimates 
deduced from the position of the lensing arc. 

distribution and our knowledge of the galaxy orbits are 
limited for studies of cluster dynamics, which allows a 
wide range of possible parametric functions for the spatial 
galaxy density distribution without even attempting the 
non-parametric methods of Merritt and Tremblay (1994) 
or considering any anisotropic orbits. This can be im- 
proved by a deep wide-field observation, to extend the 
galaxy number density distribution to a large radius (up 
to 3 Mpc) and to allow a direct measure of the cluster mass 
from a weak shear analysis. This would allow us to defini- 
tively address the issue of whether or not the cluster is in 
dynamical equilibrium and constrain the range of possi- 
ble total mass distributions allowed by the wide-range of 
data from lensing effects to X-rays. In order not to bias 
the results and incorporate a wide-range of the possible to- 
tal mass density distributions, a non-parametric method 
should also be employed. With the launch of XMM and 
Chandra, we will soon able to obtain a temperature profile 
and probe the X-ray emission at the edge of the cluster 
which is crucial to the improvement of the X-ray mass 
estimates. 
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Table 2. Spectral analysis of Abell 2104 
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-03:20:46.7 





2004(*) 


153 


2 


18.83 


0.01 


20.71 


0.01 


N 


321 


15:40:12.0 


-03:20:21.1 





0706 


213 


2 


19.33 


0.02 


21.21 


0.03 


N 


323 


15:40:15.0 


-03:16:48.0 





1535 


190 


2 


19.93 


0.02 


22.74 


0.06 


Y 


324 


15:40:16.6 


-03:19:45.8 





1635 


154 


1 


19.97 


0.02 


22.36 


0.05 


Y 


325 


15:40:17.2 


-03:21:00.7 





1531 


177 


2 


17.69 


0.01 


20.42 


0.02 


Y 


327 


15:40:19.4 


-03:20:42.4 





1503 


173 


1 


18.10 


0.01 


20.71 


0.03 


Y 


328 


15:40:20.8 


-03:18:15.1 





2849 


181 


2 


19.38 


0.01 


22.33 


0.06 


N 


329 


15:40:22.6 


-03:18:14.8 





1557 


251 


1 


20.82 


0.03 


23.01 


0.07 


Y 



Notes: 

Column 1: internal reference number to Fig. ^. 

Column 2 & 3: RA and Dec (J2000). Galaxy positions were determined from the R image and should have an accuracy of ~ 
0'.'7 rms. 

Column 4: redshift 

(*) signifies the presence of emission lines: 

102: H/3, Ha, [Nil], [Sll] 

117: H/3, [Om], Ha, [Nil], [Sll] 

214: H/3, Ha 

313: H/3, Ha, [Nil], [Sll] 

320: He I, Ha, [Nil], [S n] 

214: [On], [Om], H/3 , Ha, [Nil], [Sn] 
Column 5: AV is the internal measurement error and is related to the correlation coefficient C COT r by the formula 
AV = k/(l + Ccorr) where k ~ 870 kms" 1 was determined by the Tonry and Davis (1979) method. 
Column 6: redshift measurement quality: 

1: highest peak in the correlation function and checked by hand, 

2: highest peak in the correlation function but unable to be checked by hand. 
Column 7, 8, 9 & 10: R, AR, B and AB magnitudes: 

Column 11: Cluster member galaxy (within ±3000 km/s of the cD galaxy). 



